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Abstract Debris disks are optically thin, almost gas-free dusty disks observed around 
a significant fraction of main-sequence stars older than about 10 Myr Since the circum- 
stellar dust is short-lived, the very existence of these disks is considered as evidence that 
dust-producing planetesimals are still present in mature systems, in which planets have 
formed — or failed to form — a long time ago. It is inferred that these planetesimals 
orbit their host stars at asteroid to Kuiper-belt distances and continually supply fresh dust 
through mutual collisions. This review outlines observational techniques and results on 
debris disks, summarizes their essential physics and theoretical models, and then places 
them into the general context of planetary systems, uncovering interrelations between the 
disks, dust parent bodies, and planets. It will be shown that debris disks can serve as trac- 
ers of planetesimals and planets and shed light on the planetesimal and planet formation 
processes that operated in these systems in the past. 
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1 INTRODUCTION 

An inventory of our own planetary system uncovers its complex architecture. Eight known planets are 
arranged in two groups, four terrestrial ones and four giants. The main asteroid belt between two groups 
of planets, terrestrial and giant ones, comprises pla netesimals that failed to grow to p lanets because of 
the strong perturbations by the nearby Jupiter (e.g., ISafronovi \l96^ IWetherillL ll98oi) . The Edgeworth- 
Kuiper Belt (EKB) exterior to the Neptune orbit is built up by plan etesimals t hat did not form plan- 
ets because the density of the outer solar nebula was too low (e.g., ISafronovi Il969t iLissaue E [19821 
iKenvon & Bromlevl l2008h . Both the asteroid belt and the Kuiper belt are heavily sculptured by plan- 
ets, predominantly by Jupiter and Neptune respectively. Accordingly, they include non-resonant and 
resonant families, as well as various objects in transient orbits ranging from detached and scattered 
Kuiper-belt objects through Centaurs to Sun-grazers. Short-period comets, another tangible population 
of small b odies in the i nner s olar system, must be genetically related to the Kuiper belt that act as their 
reservoir (lOuinn et al.lll99()l) . Asteroids and short-period comets together are sources of interplanetary 
dust, obse rved in the planeta ry region, although their relative contribution to the dust production remains 
uncertain dOriin et allbOOlh . And this complex system structure was likely quite different in the past. It 
is argued that the giant planets and the Kui per belt have originally formed in a more compact configura- 
tion (the "Nice model". iGomes et al ]l2005h . and that it went through a short-lasting period of dynamical 
instability, likely explaining the geologically recorded event of the Late Heavy Bombardment (LHB). 

Similar to the s olar system, p lanetary systems around other stars are more than a star itself and one or 
several planets. As lWvattI (|2008|) pointed out, at the end of the protoplanetary phase a star is expected to 
be surrounded by one or all of the following components: various planets from sub-Earth to super- Jupiter 
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size; remnant of the protoplanetary disk, both dust and gas; planetesimal behs in which solids continue 
to grow; planetesimal belts, which are being ground down to dust. And indeed, many systems comprise 
numerous smaller objects, ranging in size presumably from dwarf planets (like Pluto or Ceres) down 
to dust. Evidence for this comes from observations of the thermal emission and stellar light scattered 
by that dust. A common umbrella term for all these "sub-planetary" solids is a "debris disk". Debris 
disks are aftermath of planet formation in the past, and they formed in the same process as the planets 
did. Even in mature systems where the planet formation has long been completed, they continue to 
evolve collisionally and dynamically, are gravitationally sculptured by planets, and the dust they produce 
through ongoing collisional cascades responds sensitively to electromagnetic and corpuscular radiation 
of the central star Hence debris disks can serve as indicators of directly invisible small bodies, planets, 
and even stars, are tracers of their formation and evolution, and represent an important component of 
planetary systems. 

We start with outlining observational techniques and results on debris disks (Section |2l) and sum- 
marizing their essential physics (Section |3]l and models (Section |4]i. In the subsequent sections, the 
observational data and their theoretical interpretation are used to draw conclusions about the "dust end" 
of the size distribution, dust itself (Section [Sj, its parent planetesimals (Section |6l), planets expected to 
be present in debris disks systems (Section |7]i, and the entire planetary systems (Section [8]l. Section |9] 
provides a summary and Usts open questions of the debris disk research. 



2 OBSERVATIONAL METHODS 



2.1 Infrared Excesses 



An efficient, and historically the first successful, way of detecting circumstellar dust is the infrared 
(IR) photometry. If dust is present around a star, it comes to a thermal equilibrium with the stellar 
radiation. With equilibrium temperatures of dust orbiting a solar-type star at several tens of AU be- 
ing typically several tens of Kelvin, dust re-emits the absorbed stellar light at wavelengths of sev- 
eral tens of micrometers, i.e. in the far-IR. As a result, the dust far-IR emission flux may exceed the 
stellar photospheric flux at the same wavelength by two-three orders of magnitude. The first discov- 
ery of the excess IR emission around a main-seque nce star, namely Vega, was done by the Infrared 
Astronomical Satellite (IRAS) (lAum ann et al 1 IT981 . Since then, IR surveys with IRAS, ISO, Spitzer, 
and other space-based and ground-based telescopes have shown "the Vega phenomenon" to be common 
for main-sequence stars. In particular, various photometric surveys of nearby sta rs have been conducted 
with the Spitzer Space Telescope. These are the GTO survey of FGK sta rs (iBeichman et al.l 2005 



Brvden et al., 2006; Bei chmag et al., 2006), th e FEPS Legacy project (.Mever et al.1, 120041: 

sll 



2005 



Hillenbrand et al.l boOSi: ICarpenter et"al.i |2009|) flie A star GTO program s jRieke et al.l 12005 



Su et al 



Kim et al 



2006h . youn g cluster programs jGorlova et al.l 120061; I Siegler et aUl2007h . binary star programs 



(iTrilling et al.ll20(j7h . as well as their combinations (e.g. lTrilling et al.U2008h . These observations were 
done mostly at 24 and 70 fim with the MIPS photometer and resulted in hundreds of detections. At 
MIPS sensitivity level, the incidences of debris disks around different-ages stars w i th spectral classes 
from A to K is about 15% (ISu et al.1 l2006l ISiegler et"an. l2007l iTYilUng et al.1 l2008t iHillenbrand et all 
l2008b . The PACS instrument aboard the Herschel Space Observatory is now performing observa- 
tions (photometry and low-resolution "small maps") at a much higher sensitivity level than that of 
Spitzer/MIPS. The sub-mm wavelength region is being probed by CSO/SHARC, Herschel/SPIRE, 
JCMT/SCUBA, APEX/LA BOCA, and o ther instruments. Soon more high-quality data should start com- 
ing from JCMT/SCUBA-2 (iMatthews et al...2007.) . 



2.2 Imaging 



So far, about 20 brightest debris disks have been resolved at various wavelengths from visual through 
mid- and far-IR to sub-mm and radio; see, e.g., http://www.circumstellardisks.org Scattered light im- 
ages in the visual and near-IR are being taken, for instance, with HST/ASC and many ground-based in- 
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struments in the coronagraphic mode, thermal emission images in the mid-IR with Gemini South/TReCS 
and Keckll/MIRLIN, whereas the JCMT/SCUBA instrument has been particularly successful in re- 
solving debris disks in the sub-mm. Resolved systems span ages from about 12 M yr for the members 
of the Beta Pictoris M oving Group dZuckerman et al.l BOOll lOrtega et al. , 20021) such as Pic itself 
( [Smith & Terrilell 1984 or AU Mic (iKalas et al.Ll2004 to - 4-12 Gvr (iDi Folco et allllool for r Get 
(iGreaves et al. r i2004 . All of the disks resolved so far exhibit various structures: rings, inner voids, 
clumps, spirals, warps, wing asymmetries. An interpretation of these structures is dicsussed in Sect.|7J 



2.3 Infrared Spectroscopy 

Early successful observations of several debris disks were made with NTT/TIMMI2 dSchiitz et al.L 
l2005b . However, most of the information comes from the Spitzer/IRS instrumen t that have taken g ood- 
quaUty spectra of more than a hundred of stars from 5 to 35 fim. For instance, iGhen et all ( |2006|) de- 
scribe IRS spectra of 59 stars with IRAS 60 /im IR excesses, five of which show a clear silicate feature 
at 10 /im. Often other spec tral features are se en that may be indicative of Fe-rich sulfides, water ice, 
and other compounds (e.g. iLisse et al.l l2007l) . In addition, another Spitzer instrument, MIPS, offered 
the so-called spectral energy distribut ion (SE P) mode, in which low-resolution spectra in the region 
from 55 /im to 90 /im have been taken dChen et al.n2008i) . In the future, high-quality mid-IR spectra are 
expected, e.g., from JWST/MIRI. 



2.4 Other Methods 



One method is near-IR polarimetry. A certain degree of linear polarization of scattere d light coming 
from circumstellar debris dust is expected, because it has a flattened, disk geometry (e.g. lKrivova et al.i 
I2OOOI) . S uccessful po larimetry observat ions were done s o far only for a few brightest debri s disks, such 
as /3 Pic dGledhill et al., 1991; Wolstencroft et al.lll995h and AU Mic (iGraham et al.ll2007h . Such mea- 
sureme nts can be particularly useful to further constrain the size distribution and chemical composition 
of dust (lKrivovaetal.Ll2000l) . Another method is near-IR interferometry, which probes dust in extreme 
proximity to the stars (the so-called "exozodiacal clouds"). Using this method, a handful of exozodis 
within 1 AU from the stars was discovered in the recent years with the GHARA/FLUOR instrument in 
the K band at 2.2 /tm (lAbsil et all I2OO6I iDi Folco et al.l l2007l: lAbsil et al.l 120081: lAkeson et al.l l2009b 
and with the Keck Interferometer Nuller at 10 /im (IStark et al.ll2009l) An ex ample of more exotic meth- 
ods suggested to detect debris disks is microlensing ( Heng & Keetonll2009h . 



3 BASIC PHYSICS 



Every debris disk is composed of solids, spanning a huge size range from hundreds of kilometers (large 
planetesimals) down to a fraction of a micrometer (fine dust). These objects orbit the central star and 
occasionally collide with each other We start with a brief characterization of their individual orbits and 
then move on to the mutual collisions. 



3.1 Stellar Gravity and Radiation Pressure 

The force keeping a planetesimal or a dust grain on a closed orbit is the central star's gravity: 

where G is the gravitational constant, A/* the stellar mass, r the radius vector of the object and 771 
its mass. For the time being, we neglect further possible gravitational forces by planets in the system 
(see Sect. |2]for a discussion of their effects). At dust sizes (s < 1 mm), the solids are also subject to 
radiation pressure caused by the central star Like stellar gravity, radiation pressure scales as a reciprocal 
of the square distance to the star, but is directed outward. Thus the two forces can be combined into a 
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"photogravitational" force (iBurns et all 1 19791) that can be written as 



(2) 



Here, (3 is the radiation pressure to gravity ratio that depends on the grain size and optical properties. 
The simplest assumption is that of compact grains of spherical shape that are characterized by their 
radius s, bulk density p, and radiation pressure efficiency Qpi . The last coefficient controls the fraction 
of the momentum transferred from the infalling radiation to the grain. Its values range from for perfect 
transmitters to 2 for perfect backscat terers. An ideal abs orber's radiation pressure efficiency equals unity. 
In that case, the /3-ratio is given by dBurns et al.Lll979h 



/3 = 0.574 ( ^ 



Me 



1 gem 



1 /im 



(3) 



where L/Lq and AI /Mq are luminosity and mass of the star in solar units. 

If, by fragmentation or any other erosive process, smaller particles are released from larger ones, 
they start to "feel" a photogravitational force rather than gravitational force. The smaller the grains, the 
more the radiation pressure they experience compensates the central star's gravity. Thus the orbits of 
small particles differ from those of the parent bodies (Fig.[T]left). A parent body on a circular orbit, for 
example, releases fragments into bound elliptic orbits with larger semimajor axes and eccentricities up 
to /3 ~ 0.5. If 0.5 < /3 < 1, the grain orbits are hyperbolic and thus unbound. For parent bodies in 
elliptic orbits, the boundaries between bound and unbound orbits are somewhat smeared, because the 
type of orbit depends on the ejection point. All grains with /3 < 1 orbit the star on Keplerian trajectories 
at velocities reduced by a factor of \/l — j3 compared to macroscopic bodies, which are purely under 
the influence of gravitation. Finally, below the critical size where /3 = 1 the effective force is repelling, 
and the grains move on "a nomalous"(i.e. , open outward) hyperbolic orbits. Borrowing the terminology 
from solar system studies (IZook&BergLll975h . we may call dust grains in bound and unbound orbits 
a-meteoroids and /3-meteoroids, respectively. 



Hyperbolas 



Anomalous 
hyperbolas 




|3>0.10 |3>0-20 |3>0.25 




Fig. 1 Left: Three possible types of orbits of dust particles under the combined action of 
stellar gravity and direct radiation pressure. For illustrative purposes, grains are assumed to 
be released from a circular orbit. Right: Schematic of a debris disk produced by a parent 
planetesimal belt. 
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3.2 Collisions 



As the very name "debris disks" suggests, they should incessantly produce "collisional debris", so 
that destructive collisions must be a dominant process operating in these systems. For collisions to 
be destructive — and even to occur at sufficient rates - a certain minimum level of relative velocities 
is necessary. In protoplanetary disks, this condition is not fulfilled, as relative velocities are strongly 
damped by dense ambient gas, so that the orbits all have low eccentricities and inclinations. In con- 
trast, debris disks are gas-poor, and damping is not efficient at all. However, the absence of damping 
is necessary, but not sufficient for relative velocities to be high. As a protoplanetary disk transforms to 
a debris disk in the course of its evolution, solids are expected to preserve a low velocity dispersion 
they had before the primordial gas re moval. Accord ingly, to ignite destructive collisions, debris disks 
must be "stirred" by some m echanism (IWvattll2008h . This can be self-stirring by largest planetesimals 



^ , ^ Pi 

( Kenvon & Bromlevi 2004ah or stirring by planets orbiting in the inner gap of the disk dWvatlL I200I 



iMustill & WvatT 120091) . Further possibi ities which, however, are not considered typical for the ma- 
jority of the disks, include "fast ignition" by the abr upt injection of a plan et (suggestively of Neptune 
masses) into the disk after planet-planet scattering dRavmond et al.l l2009l) or external events such as 
stellar flybys dKenvon & Bromlevi I2OO2I) . 

Once the disk is sufficiently stirred, the collisional cascade sets in. The material is ground all the 
way down to the dust sizes, until the smallest fragments with p 0.5 are blown away from the disk 
by radiation pressure. Note that at dust sizes, stirri ng is no longer req uired. As typical eccentricities 
of grains are of the same order as their /3-ratio (e.g. lBurns et al.l I979I) . radiation pressure ensures the 



impact velocities to be sufficiently high. What is more, some of the hyperbolic grains, /3-meteoroids, 
may be energetic enough to tr i gger a self-sustained casc ade ("dust avalanche") on their way outward 



through the disk dArtvmowiczl I997tlKrivov et al" , 2000l). This p rocess, however, can only be efficient 



in the dustiest debris disks such as j3 Pic dGrigorieva et 



y) inispr 
dl l2007f) . 



We now take a more quantitative look at the collisions. An important quantity in the collisional 
prescription is the critical specific energy for disruption and dispersal, Qq. It is defined as the im- 
pact energy per unit target mass that results in the largest remnant containing a half of the original 
target mass. For small objects, Qq is determined solely by the material strength, while for objects 



scribed by the sum of two power laws (see 


e.g., 
19981; 


Davis et all 19851 Holsappld. 1994; Paolicchi et al. 


19961 Durda & Dermottl 1997HDurda et al.. 


Benz & AsDhaud.ll999HKenvon & Bromlevi l2004b 


Stewart & Leinhardtll2009l): 



Qt) — Qs 



s 
1 m 



1 km 



(4) 



where the subscripts s and g stand for the strength and gravity regime, respectively. The values of 
Qs and Qg lie in the ra nge ^ 10^-10^ erg g~^, 6s is between and 0.5, and bg between 1 and 2 
(iBenz & Asphauglll999l) . The critical energy reaches a minimum (^ 10^-10^ erg g^^) at sub-km sizes. 
At dust sizes, Eq. dU suggests 10® ergg~^, but true values remain unknown, as this size regime has 
never be en probed experime ntally. Both colliders are disrupted if their relative velocity Uroi exceeds 
(see, e.g. lKrivov et al.ll2005l their Eq. 5.2) 



/2(mt +mp)^ 



mtmp 



(5) 



where mt and nip are masses of the two colliders. Thus two like-sized objects will be destroyed if they 
collide at a speed exceeding y^SQ^. For two dust grains, the critical speed is several hundreds m s^^. 
At several tens of AU from a star, this would imply eccentricities of ~ 0.1 or higher Such eccentricities 
can be naturally attained by dust grains with /3 ~ 0. 1 as a result of the radiation pressure effect described 
above. 
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The actual physics of colUsions in debris disks is more compHcated. ColHsions occur across the 
size range that is ^ 12 orders of magnitude wide and at quite different velocities. Statistically, collisions 
between objects of very dissimilar sizes (e.g., impacts of dust grains onto "boulders") outnumber those 
between like-sized objects. Besides, grazing or at least oblique collisions are more frequent than frontal 
ones. All this leads to a spectrum of possible collisional outcomes that take place simultaneously in the 
same disk, comprising further growth of large planetesimals, partial destruction (cratering) of one or 
both projectiles with or without reaccumulation of fragments, as well as a complete disruption of small 
planetesimals and dust particles. Detailed analyses sho w that two outcomes are expected to play a crucia l 
role in debris disks: complete disruption and cratering dThebault & Augereaull2007tlMuller et alll2010l) . 



3.3 Drag Forces 

A more accurate, special-relativistic, derivation of the photogravitational for ce leads to an add itional 
velocity-dependent term that has to be added in the right-hand side of Eq. ^ dBurns et al.Lll979l) : 



-r PR — o 



/vr\ r V 
— - + - 

\cr J r c 



(6) 



with V being the velocity vector of the particle, which is referred to as the Poynting-Robertson (P-R) 
force. Being dissipative, this force causes a particle to lose gradually its orbital energy and angular 
momentum. Th us all g-meteoroids mov e in Keplerian ellipses with reduci ng semimajor axes a an d 
eccentricities e (IWvatt&Whipplelll950h (except very close to the central star. lBreiter & Jacksonll998h . 
On timescales of thousands of years or more the trajectory shrinks to the star 

Like stellar electromagnetic radiation gives rise to radiation pressure forces, stellar particulate radi- 
ation — stellar wind — causes stellar wind forces. Similarly to the net radiation pressure force, the total 
stellar wind force can be decomposed to direct stellar wind pressure and stellar wind drag. For most of 
the stars, the momentum and energy flux carried by the stellar wind is by several orders of magnitude 
smaller than that carried by stellar photons, so that the direct stellar wind pressure is negligibly small. 
However, the stellar wind drag forces cannot be ignored, because the stellar wind velocity Vgw, to replace 
c in Eq. (|6]l, is much smaller than c. Stellar wind force s can be very important in debri s disks around 
late-type stars jPlavchan et al.L 120051: [Strubbe & Chiangll2006l: lAugereau & Beuslll2006h . 

Compared to gas-rich protoplanetary disks, debris d isks are expected to be extremely gas-poor, be- 
cause primordial gas disperses on timescales ^ 10 Myr d Alexander, ,2008; .Hillenbrand. .2008,) . Even so, 
in youngest debris disks with ages of 10-30 Myr, little amounts of gas have been detected. Debris 
disk gas might be either the remnant of the primordial protoplanetary disk or, more likely, of sec- 
ondary origi n. Possible mechanisms for producing secon dary gas include photon-ind uced desorption 
from solids (IChen et al.l |2007[) and grain-grain collisions dCzechowski & Mannll2007l) . Part of the ob- 
served gas may als o stem from comet evaporation, as inferred from observed time-variable stellar ab- 
sorption lines (e.g. Ferlet et allll987t lBeujt_&_yaliTOn. 2007). In some systems, exemplified by /? Pic 



(lOlofsson et all 1200 ih and HP 32297 dRedfieldl I2007L evidence for a stable, orbiting gas component 



was found. Dy namical effects exerted by rot ating gas on dust in youn g debris disks have been studied, 
for instance, bv lThebault & Augereaul (l2005h and lKrivov et"ai1d2009l) . 



3.4 Other Forces and Effects 



Dust grains interacting with the stellar radiation and stellar wind acquire electrostatic charges due to a 
variety of effects (photoelectron emission, secondary electron emission, sticking of electrons and ions 
etc.). Under the presence of the star's magnetic field, they experience the Lorentz force. Its influence on 
the dynamics of interplanetary dust particles in the solar system is well understood. As the dust grains 
move through the sectored magnetic field of the Sun with alternative polarities, the Lorentz force rapidly 
changes its direction. The mutual near-cancellation of these contributions is not complete, however On 
long time spans the Lorentz force result s in a "Lorentz diffusion", i.e. essentially stochastic changes in 
a, e and, most importantly, inclination i dMorfifl & Grunill979l:IConsolmagnolll980tlBarge et"an,ll982l: 
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IWallis&Hassan[ll985h . The effect becomes significant for grains less than several micrometers in size. 
The Lorentz force is not included in most of the debris disk models. This is because it affects relatively 
small grains, which make only a minor contribution both into the mass budget and the total cross section 
area of the dust disk. 

Sublimation, or transition of dust grains from the solid to the vapor state, occurs in the vicinity of the 
star. The radius of the sublimation zone depends strongly on the material and porosity of the dust grains 
and on the luminosity of the prim ary. In the solar system, silicate and carb onaceous grains sublimate at 
2-4 solar radii from the Sun (e.g. lKimura et al.lll997tlKrivov et al. I [19981) . However, for a water ice - 
organics mixture the sublimation distance in the solar system is as large as ^ 20 AU, and increases to 
30-40 AU around luminous A-stars (Kobayashi et al.. 20()9i)- 

A nu mber of more exotic effects exist that may influence the dust particles themselves and their 
dynamics dBurns et al.Lll979l : lKapisinskvlll984l) . Some are related to the fact that the star is not a point- 
like object (differential Doppler effect, radiation pressure from an extended source). Others stem from 
stellar radiation together with rotation of objects (Yarkovsky effect, windmill effect, Radzievsky effect). 
Still others are related to modification of the grain properties by the environment (sputtering by plasma 
or interstellar grains, packing effect, electrostatic breakup), etc. None of these, however, seem to be 
of generic importance for debris disks, requiring either closeness to the star, smallness of particles, or 
special material compositions. 



4 MODELS 



4.1 Numerical Simulation Methods 

As described in Sect. [3] each debris disk can be treated an ensemble of objects of different sizes (from 
dust to planetesimals) that orbit the star under the influence of gravity, radiation pressure, and other 
forces. In addition, they experience collisions which destroy or erode these objects and produce new 
ones, collisional fr agments. Such a system can be numerically modeled by a variety of methods (see 
iKrivov et al 1 12001 for a summary) that can be classified into three major groups: iV-body simulations, 
statistical approach, and hybrid methods. 

N-body codes follow trajectories of individual of disk objects by numerically integrating their equa- 
tions of motion. During numerical integrations, instantaneous positions and velocities of particles are 
stored. Assuming that the objects are produced and lost at constant rates, this allows one to calculate 
a steady-state spatial distribution of particles in the disk. Then, collisional velocities and rates can be 
computed and collisions can be applied. It is usually assumed that each pair of objects that come in 
contact at a sufficiently high relative velocity is eliminated from the system without gener ating smaller 
fragments dLecavelier des Etangs et all 11996! : Stark & Ku chneii 120091 : iKrivov et al.L 20091) or produces 



a certain number of fragments of equal size (iWyatti i2006l) . The TV-body method is able to handle an 
arbitrary large array of forces and an arbitrary complex dynamical behavior of disk solids driven by 
these forces. Accordingly, it is superior to the others in studies of structures in debris disks arising from 
interactions with planets, ambient gas, or interstellar medium. However, this method cannot treat a large 
number of objects sufficient to cover a broad range of particle masses and thus is less successful in 
modeling the collisional cascade. In particular, an accurate characterization of the size distribution with 
A^-body codes is hardly possible. 

Statistical method effectively replaces particles themselves with their distribution in an appropri- 
ate phase space. The ideas trace back to the classical calculation of the velocity distribution in par- 
ticle ensembles with Boltzmann equation (Boltzmann, 1896; Chapman & Co wling. 1970) and of th e 
mass (or size) distribution with the coagulation equation dSmoluchowskii 1 1 9 1 61; Ichandrasekhail 1 1 943h . 
Note that the term "coagulation equation" is actually used regardless of whether the colliding parti- 
cles merge, fragment, or get destroyed. In applications to planetesimal and debris dis ks, one intro 



duces a mesh of several variables comprising, for i nstance, mas s, distance, and velocity dSpaute et al. 
199ll:IWeidenschining et allll997HKenvon & Luul [l99 8[ll999allbH K" envon & Bromlevll2002i l2004allc : 



Thebault et al.ll2()()3tlThebault & Augereaull2007l) or mass and orbital elements dPelFOro & Paolicchi , 
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1 19981 iDeir Pro et al 1 1 19981 1200 ll l2002t iKrivov et all l2005h . The number of particles in each bin of 
the mesh at successive time instants is computed by solving equations that describe gain and loss of 
objects by collisions and other physical processes. Many of the results presented in this paper (e.g.. 
Fig. |2|i were obtained with the code ACE (Analysis of Collisional Evolution) that simul ates the disks 
in such a manner, using a 3D mesh of masses, pericentric distances, and eccentricities ( iKrivov et all 
I2005ll2006ll2008l) . Statistical codes are much more accurate in handling collisions than iV-body ones, 
but treat dynamics in a simplified way, for instance, by averaging over angular orbital elements. For this 

reason, they are less suita ble for simulation of structures in debris disks. 

Lastly, hybrid codes dWeidenschilling et all Il997t iGoldreich et all |2004|) combine iV-body inte- 
grations of a few large bodies (planets, planetary embryos, biggest planetesimals) with a statistical sim- 
ulation of numerous small planetesimals and dust. Su ch codes, originally developed for planetesimal 
accretion, have also been applied to debris disks (e.g. iBromlev & Kenvonl I2OO6I : iKenvon & Bromlevl 
I2OO8). Finall y, a stand-alone "superp article" method developed to study major collisional break-ups in 
debris disks (iGrigorieva et all |2007|) should be mentioned. This method, too, combines TV-body and 
statistical approaches, although in a different way. 

4.2 A "Standard Model" for a Debris Disk 

A standard model of a debris disk can be devised by considering two major effe cts, stellar pho- 
togravity and collisions, and neglecting the drag forces and all other processes (e.g. Thebault et al. , 
2003; Krivov et al 20()6; _Slrubbe_& Chiang, 2006; Thebault & Augereau, 2007; Lohne et al., 2008; 
Krivov et all I2OO8I; iThebault & Wui 20081) . Imagine a relatively narrow belt of planetesimals ("birth 
ring" or "parent ring") in orbits with moderate eccentricities and inclinations, exemplified by the clas- 
sical Kuiper belt in the solar system (Fig. [1] right). The planetesimals orbiting in the birth ring undergo 
collisional cascade that grinds the solids down to dust. At smallest dust sizes, stellar radiation pressure 
effectively reduces the mass of the central star and quickly (on the dynamical timescale) sends the grains 
into more eccentric orbits, with their pericenters still residing within the birth ring while the apocenters 
are located outside the ring. As a result, the dust disk extends outward from the planetesimal belt. The 
smaller the grains, the more extended their "partial" disk. The tiniest dust grains, for which the radiation 
pressure effectively reduces the physical mass by half, are blown out of the system in hyperbolic orbits. 
The radiation pressure blowout of the smallest collisional debris represents the main mass loss channel 
in such a disk. 

One expects a balance between the production of dust by the collisional cascade and its losses by 
radiation pressure blowout. If such a balance exists, the amounts of particles with different sizes on 
different orbits stay constant relative to each other, and a debris disk is said to be in a quasi-steady 
state. However, the absolute amounts should adiabatically decrease with time, because the material at 
the top-end of the size distribution is not replenished (therefore "quasi"). For brevity, "quasi" is often 
omit ted and a simple "steady state" is used. Basic properties of steady-state disks, derived analytically 
(e.g. Strubbe & Chiang, 2006) and by numerical simulations with statistical codes dKrivov et al.l l2006t 
rfhebault & Augereau., ,2007.) are described in subsequent subsections. The steady-state evolutionary 
regime can be perturbed, for instance, by occasional collisional break-ups of largest planetesimals or by 
"shake-downs" of the system due to instability of nearby planets. After such events, the disk needs some 
time to relax to a new steady state. 

4.3 Collision-Dominated and Transport-Dominated Disks 

The model just described is only valid, if the collisional timescale is shorter than the characteristic 
timescale of the drag forces. Such disks are usually referred to as collision-dom inated, as opposed to 
systems that — at dust sizes — are transport-dominated (e.g. IKrivov et al I I2OOOI) . In the latter case, ra- 
dial transport of dust material by various drag forces occurs on shorter timescales then collisions. Then, 
additional removal mechanisms may play a significant role. For example, Poynting-Robertson (P-R) 
drag can bring grains close to the star where they would sublimate or deUver them into the planetary 
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region where they would be scattered by planets. To be transport-d ominated, the system should either 
have an optical depth below current detection limits (IWvatil2005b or be subje ct to transport mecha- 
nisms other than P-R drag, such as strong stellar winds around late-type stars dPlavchan et all 120051 : 
IStrubbe & ChiangibOOetfAugereau & B"eus3.l2006l) . Most of debris disks detected so far are thought to 
be collision-dominated. Accordingly, transport-dominated systems will not be considered here. 

Note that at sufficiently large sizes, all disks are dominated by collisions. This is because the lifetime 
against catastrophic collis ions in a disk with the s~'^-^ size distribution discussed below scales as cx -y/s 
(see, e.g.. Appendix A in lWvatt et al.l [l999l) . whereas any drag force is proportional to the ratio of the 
mass and cross section, i.e. oc s. For instance, the interplanetary dust cloud in the solar system at 1 AU 
from the Su n is tenuous enough to be P-R-dominated at dust sizes, but is collision-dominated above 
s « 100 /im jGriin et all [19851) . 



4.4 Size and Radial Distribution in Collision-Dominated Disks 

We now consider typical size and radial distribution of material in a steady-state collision-dominated 
debris dis k, seen in simulation s. As an example, we take the "reference model" of the debris disk of 
Vega from lMulleretal.1 (1201 Ol) that assumes a parent ring between sa 80 and « 120 AU. However, the 
results are generic and qualitatively hold for all collision-dominated disks. 



Size distribution „ Radial distribution 




u 1^m Imm 1m 1l<nn u 60 80 100 150 200 300 500 

radius distance r[AU] 



Fig. 2 Modeled dust distribu tions in a collision-dominated debris disk, here: the "reference" 
model of iMuller et al] ( |2010|) . Left: size distribution at different distances from the star with 
catastrophic and cratering collisions included in the model (thick lines). For comparison, a 
distribution at 1 00 AU computed without cratering is also shown (thin solid line). Straight 
line depicts the iDohnanyil /l969l) distribution for comparison. Right: radial distribution of 
different-sized grains (thin lines), as well as the overall radial profile of the normal geometrical 
optica l depth (think solid line). Straight line is the r^^-^ analytic prediction for the optical 
depth dStrubbe & Chiangll2006h . 

Figure |2] (left) depicts the size distribution. Most noticeable is the peak at the particle size that 
corresponds to /? ~ 0.5, where the radiation pressure is half as strong as gravity. For the Vega disk 
this blowout radius is Sbiow ~ 4 /im. Below this size bound orbits are rare (and below the size where 
/3 = 1 they are impossible), and the grains are blown away on a disk-crossing timescale of the order 
of 10^ to 10'^ yr. Typically, the amount of blowout grains instantaneously present in the steady-state 
system is much less than the amount of slightly larger grains in loosely bound orbits around the star 
This is because the production rate of the grains of sizes slightly less and slightly greater than Sbiow 
is comparable, but the lifetime of bound grains (due to collisions) is much longer than the lifetime 
of blowout grains (disk-crossing timescale). This explains a drop in the size distribution around the 
blowout size seen in the left panel of Fig. |2] The size distribution above Sbiow sensitively depends 
on the collisional outcome model assumed, which is seen by comparing the distributions at 100 AU 
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computed with and without taking cratering colUsions into account. The overall size distribution inside 
the parent ring (the 100 AU curve) at sizes s ^ Sbiow roughly follows the classical Dohnanyi power law, 
which cor r espon ds to the differential size distribution n cx s~^'^. That law was theoretically derived by 
lDohnanvil(ll969l) for a collisional cascade in a system without sources and sinks, composed of objects 
with = const. The deviations of the size distribution in Fig. |2] (left) from the Dohnanyi law reflect 
the size dependence of Q^, given by Eq. (|4]i. A notch at s ~' 100 m corresponds to the minimum of 
Qq, and the slope at larger sizes is slightly s teeper that Dohnanyi's, because increases with size in 
the gravity regime jPurda & DermotA Il997l) . Another generic feature of the size distribution is that it 
becomes narrower at larger distances from the star. At twice the distance of the parent ring, 200 AU, the 
distribution converts to narrow peak around ~ 2sbiow composed only of small, high-/?, barely bound 
grains sent by radiation pressure into eccentric orbits with large apocentric distances. 

The spatial distribution of dust in the same modeled disk is illustrated by Figure |2] (right). As noted 
above, smaller grains with higher (3 ratios acquire higher orbital eccentricities. Since the eccentricities 
of particles slightly above the blowout limit (s = 5 fim) are the highest, their radial distribution is the 
broadest, whereas larger particles stay more confined to their birth regions. In the figure, this effect can 
be seen from how the curves gradually change from the smallest (s = 5 /xm) to the largest bound grains 
(50 /im). The latter essentially follow the distribution of the parent bodies. In addition, the same figure 
shows the radial profile of the nor mal geometrical opt i cal de pth. The slope of r oc r~" is close to 
a = 1 .5, as predicted analytically bv lStrubbe & Chiang (|2006|) for collision-dominated disks. 



4.5 Disk Appearance at Different Wavelengths 

Size segregation of dust in debris disks described above means that in resolved images at different 
wavelengths the same disk would look differently (Fig. [3]). Measurements at longer wavelengths (sub- 
mm) probe larger grains, because they are cooler (see Sect. |5]), and thus trace the parent ring. Such 
observations may also reveal clumps, if for instance there is a planet just interior to the inner rim of 
the parent ring, and parent planetesimal s and their debris are trapped in e xternal resonances, similar to 
Plutinos in 3:2 resonance with Neptune dWvatli I2OO6I : iKrivov et allbOOTi) . At shorter wavelengths (far- 
IR, mid-IR), smaller (warmer) grains are probed. Thus the same disk appears much larg er The image 
may be featureless, even if the parent ring is clumpy, bec ause strong radiation pressure dWvati l2006l) 
and non-negligible relative velocities jKrivov et all |2007|) would liberate small particles from resonant 
clumps. As a result, they would form an extended disk, as described above, regardless of whether their 
parent bodies are resonant or not. Finally, at shortest wavelengths of thermal emission, only the hottest 
closest-in grains are evident in the observations, and thus again only the parent ring is seen. 



4.6 Scaling Laws for Evolution of Collision-Dominated Disks 

Some useful laws that describe quasi-steady state evolution of colUsion-dominated disks have been 
found dWvatt et alir2007aHLohne e"tlill2008tl&ivov et alllloosl) . 

1. Dependence of evolution on initial disk mass. Consider a narrow parent ring with initial mass 
M{t = 0) = il/o at a distance r from the star with age t. Denote by F{Mo, r, t) any quantity directly 
proportional to the amount of disk material in any size regime, from dust grains to planetesimals. In 
other words, F may e qually stand for the total disk mass, the mass of dust, its total cross section, etc. 
There is a scaUng rule jLohne et al.Ll2008l) : 

F{xMo,r, t) = xF{Mq, r, xt), (7) 

valid for any factor a; > 0. This scaling is an exact property of every disk of particles, provided these 
are produced, modified and lost in binary collisions and not in any other physical processes. 

2. Dependence of evolution on dis tance. Another scaling rule is the depe ndence of the evolution 
timescale on the distance from the star dWvatt et ali[2007aHLohne et al.Ll2008h . Then 

F(Afo, xr, t) « F(A/o, r, x-^'H). (8) 
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Fig. 3 Appearance of a debris disk at different wavelengths from 10 to 850 nva. Top row : 
synthetic images of the Vega disk, based on the "reference" model of iMiiller et alj (|2010|) . 
Bottom row: the same images convolved with the Point-Spread Function (PSF), which we 
assumed to be a Gaussian with the width {a) of 0.089" at 10 /im (Keck, mirror diameter 
10 m), 2.37" at 24 /im and 1.77" at 70 /im (Spitzer/MIPS, 0.9 m), 2.53" at 100 /im and 4.05" 
at 160 /tm (Herschel/PACS, 3.5 m), 2.66" at 450 /tm and 5.03" at 850 /tm (JCMT/SCUBA, 
15 m). These convolved images show how the Vega disk would appear in observations with 
best present-day instruments operating at t he respective wavelengths. Note that in a ctual sub- 
mm images the Vega disk appears c lumpy (Holland et al.lll998l : lMarsh et alll2006l) . which is 
not seen here, because the model bv lMiiller et al. ( 2010l) is rotationally-symmetric. Courtesy: 
Sebastian Milllen 

Unlike Eq. (|7]i, this scaling is approximate. 



3. Dust mass as a function of time. Finally, the third scaling rule found in lLohne et al.l (l2008h is the 
power-law decay of the dust mass 

F{Mo,r,xt)^x-^F{Mo,r,t), (9) 

where ^ w 0.3 . . . 0.4. This scaling is also approximate and, unlike Eq. (|7]i and Eq. only applies to 
every quantity directly proportional to the amount of dust. In this context, "dust" refers to all objects 
in the strength rather than gravity regime, implying radii less than about 100 meters. The scaling is 
sufficiently accurate for disks that are much older than the collisional lifetime of these 100 m-sized 
bodies. 



5 DEBRIS DISKS: SEEING DUST 

5.1 Dust Temperature and Fractional Luminosity 

Observational data for most of the debris disks are pretty scarce and consist only of a few photometric 
points indicative of the IR excess. In many cases, these points can be fitted with a spectral energy 
distribution (SED) of a blackbody with a single temperature. Two useful quantities that can be derived 
from such observations are the dust temperature Td and the dust fractional luminosity /d, defined as a 
ratio of the bolometric luminosities of the dust and the star. 

Both Td and /d can be est imated from the wavelength where the dust emission flux peaks, A^^^x' 
and that maximum flux, F^-^^^ dWvatll l2008l) . (Henceforth, any "flux" is assumed to be specific flux per 
unit frequency interval.) If dust behaves as a blackbody, Wien's displacement law gives: 

Td = 5100K(i^). (10) 

\ max / 

Assuming, further, that also the stellar photosphere emits as a blackbody, it is straightforward to find 
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where X^^^^ and F^^-^ are the wavelength where the stellar radiation flux has a maximum and the 
maximum flux itself, respectively. Note that Eq. (fTTT i can only be used for ballpark estimates, because 
the actual stellar spectrum may deviate from a single blackbody spectrum considerably. A more accurate 
method to determine /d that does not imply any blackbody assumptions is using its definition: one 
calculates the bolometric dust luminosity from the SED, uses the bolometric luminosity of the star from 
catalogs, and takes their ratio. 



5.2 Dust Location, Dust Mass, Dust Sizes 

Keeping the assumption that dust interacts with the stellar radiation as a blackbody and making some 
assumptions ab out the disk g eometry, a number of further useful physical parameters can be derived 
from Td and /d dWvatllbOOSl) . For instance, assuming the disk to be a narrow ring between r — dr and 
r + dr, dust temperature gives the distance from the star: 



, , (12) 

whereas /d is directly related to the total cross-section area of dust a^, 

<Jd = W/d, (13) 
and to the normal geometrical optical depth t±, 

rx = (r/dr)/d. (14) 

Dust mass Md can be derived from fractional luminosity as well. However, this estimation requires 
further assumptions, namely about the dust bulk density pd and a "typicaf'dust size Sd (which is poorly 
defined, because /d-dominating and A/d -dominating sizes are different). Then, 

16 2 

Md = yTrpdSd?' Id- (15) 

More detailed analyses go beyond the assumptions of blackbody, confinement of dust to a narrow 
radial zone, and like-sized dust grains. A usual procedure (e.g.;Wolf & H illenbrand, 2003) is to assume 
one or another chemical composition of dust, as well as a radial distribution of dust from to r,nax 
and a size distribution from Smin to Smax- Both distributions are usually postulated to be power laws. 
Both the limits and slopes of these power laws are treated as free parameters, and their values are sought 
by fitting the SED and, where avaliable, the resolved images. 

This procedure is v ery efficient, have been implemented in a number of codes (e.g. 
IWolf & Hillenbrand |2005[) and is intensively used in interpretation of the debris disk observations. 
A short summary of the essential results is given below. However, it reveals two problems that need to 
be mentioned. One is that relaxing the blackbody assumption makes dust temperature dependent not 
only of the dust location, but also on the particle sizes (Fig.|4]i. Roughly speaking, smaller grains at the 
same distance are warmer than larger ones. This leads to a fundamental degeneracy between distance 
and size, as the same SED can be produced equally well by a disk of smaller grains farther out from the 
star and by a disk of larger grains closer in. One way to break this degeneracy is to invoke information 
from resolved images, if these are available, as they show directly where most of the emitting dust is 
located. Another method is to invoke theoretical constraints, notably to assume that the minimum grain 
close to the radiation pressure blowout limit Sbiow (Sect.H]). Another problem of the fitting 
approach described above is more difficult to cope with. It roots in the assumption that radial and size 
distribution are both power laws, independent of each other. As we saw, this is not what is expected 
from the debris disk physics. Instead, one expects the size distribution to strongly depend on distance 
and conversely, the radial distribution to be very dissimilar for grains of different sizes. Additionally, 
both size and radial distribution can substantially depart from single power laws (Fig.|2]). 
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Fig. 4 Calculated temperature of different-sized astrosilicate (iLaor & Drainel 1 19931) grains 
at different distances from a G2V star. Courtesy: Sebastian Miiller. 



5.3 Outer Disks, Inner disks, and Exozodiacal Clouds 



We now come to a brief overview of essential results on dust location and masses obtained with the pro- 
cedures described above. For many of the systems, the SED is consistent with a single, radially narrow 
dust disk. The typical dust temperatures are of the orders of several tens of Kelvin, which corresponds 
to the dust disks located at several tens to a hundred AU from the star We will refer to them as "outer 
disks". Such distances readily suggest that the outer disks are likely produced by the Edgeworth-Kuiper 
belt (EKB) analogs around those stars. H owever, t he typical fractiona l luminosities of ^ 10~^ — 10""^ 
are at least by two (iVitense et all 1201 Oh to four (iBooth et al.L l2009h orders of magnitude more than 
the fractional luminosity of the presumed EKB dust in the present-day solar system. The derived dust 
masses in debris disks (up to s ~ 1 mm, best estimated from sub-mm measurements), are in the range 
10-3-10° Earth mass (e.g. lSheret et"!!! |2004 Iwilliams et all [20041: iNaiita & Williamsll2005l) . 

There are also systems that cannot be explained with a single outer dust disk. Apart from, or instead 
of, the cold dust disks, some stars have dust within ~ 10 AU. The majority of such "inner disks" may 
be associated with collisionally-evolving asteroid-like belt or sustained by the cometary activity (both 
mechanisms are known to be sources of the zodiacal cloud in the solar system). Alternatively, the fact 
that many solar-type (FGK) stars show a considerable warm excess at ages < 100 Myr, which is about 
the timescale for accretion of the Earth in the solar syste m, suggests that inner dust may be debris from 
terrestrial planet formation (iKenvon & Bromlevll2004bh . The dust masses in the inner debris disks are 
estimated to be typically in the range lO-^-lO"^ Earth mass (e.g. iKrivov et dl 120081 and references 
therein). 

Yet closer-in to the central stars, within 1 AU, "exozodiacal clouds" have been de tected with near-IR 
interferometry. The m ost prominent examples are stars h osting exo zodis are Vega (Absil et all l2006h . 



T Cet and Fomalhaut (iDi Folco etali 120071: lAbsil et al.L |2008; Ake son et al. I l2009l) 7 as well as 77 Crv 



jBrvden et all l2009bl) . The archetypical exozodiacal cloud of Vega has an estimated dust mass of 
~ 10"^ Earth mass (equivalent to the mass of an asteroid about 70 km in diameter) and a fractional 
luminosity of 5 x 10^"'. 
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5.4 Frequencies of Disks around Stars with Different Spectral Types and Ages 

In recent years, detailed statistical studies have been published that show how the parame ters of outer 



and inner dust disks depend on the spectral type of the star as well on the sys t em's age (iMever et al 
20041: iBeichman et al L I200I iKim et all 120051: Rieke et al l 120051: iBrvden et all 120061: (Beichman et al. 



2006 : Su et alfUOOe l: 'Gautier et a ll l2007l ISieder et al.l 120071 iHi llenbrand et all 120081: iMever et al. 



2008[ iTrilUng et al.l 2008: Brvde n et al.l l2009ai: ICarpenter et al.l 12009: Plav chan et airi2009l) . These 



statistics are largely based on the Spitzer/MIPS measurements at 24 and 70 i-im. Note that for A-type 
stars, outer EKB-like disks are essentially probed at both wavelengths. But for G-stars, only 70 fim 
excesses correspo nd to the EKB region while 24 /im measurements detect inner dust inside ^ 10 AU 
(see, e.g.. Fig 1 in IWvatlll2008h . Still longer wavelengths may be required to detect outer disks around 
K and especially M dwarfs, because EKB-sized disks around such stars are very cold. 

A bout 33% of A- stars over all ages show excess emission at 24 /xm and/or 70 /im, indicative of outer 
disks dSu etalil2006l) . T he disk fre quency around A-stars is as high as 60% at ^ 10 Myr and declines 
to about 10% at 600 Myr (ISiegler et al ... 20071) . In addition, d isks around younger star s are brighter, with 



a possible broad peak of the 24 fiui luminosity at ^ 10 Myr dHernandez et al 



jer star s 
" l2006h . 



For sol ar-type (FO-K O) stars, the average frequency of outer disk detections at 70 fim was found to 
be w 16% dBeichman et al ., 2006; Bryden et al., ^006; T rilling et_al,, 2008). For stars cooler than Kl, 
the fraction may drop to dBeichman et al.l 20061). and low detection frequencie s for K and M 



stars have been confirmed by other analyses ( Gautier et al.l 120071 iPlavchan et al.l l2009l) . However, as 



noted above, their detection should be more efficient at wavelengths longer than 70 /J,m. Indeed, a sub- 
mm JCMT/SCUBA sur vey for M dwarfs sug gests that the proportion of debris disks is likely higher 
than found with Spitzer dLestrade et am2006l) . The frequencies of disks around FGK stars decline with 
the age more slowly that those of A-type stars. The excess strength, or fractional luminosity, exhibits a 
smooth decay with the stellar age as well. 

All authors point out a decay of the observed mid- and far-IR excesses with the system age. 
However, the dust fractional luminosity exhibits a large dispersion at any given age even within a nar- 
row range of spectral classes. We will come to the interpretation of the long-term decay in debris disks 
in Sect. |6] Apart from the time evolution of the fractional luminosity, attempt s have been made to find 
correlations between the disk radius on age. From IRAS 60 and 100 /im data. lRhee et al. I (l2007l) found 
that the radii of largest disks increase with the age, but sim ilar analyses based on 24 and 70 /im Spitzer 
detections do not reveal any significant correlation dWvatt e t al. 2007b and T. Lohne, pers. comm). 

Inner disks with a moderate fractional luminosity around A-type stars are difficult to detect because 
it requires measurements shortward of 24 fmi where the photosphere is bright. Nevertheless, they can 
be detected through measurements at 24 yum around solar-type stars. The incidence s of inner disks drop 
from 20-40% a t ^ 20 Myr to only a few percent at ^ 1 Gyr dSiegler et all |2007'; 'Me ver et all I200S 



ICarpenter et al.l l2009h . The detection fraction averaged over the ages is ~ 4% (iTrilling et al. 



Many of older 100 Myr) FGK stars with inner dust have an abnormally high fractional luminosity, 
as discussed in Sect. |6l Apart from the mid-IR photo metry, evidence for i nner dust can also be found 
in the Spitzer/IRS spectra in the « 5 to 35 /im range. lEawler et all d2009b analyzed IRS observations 
of nearby solar-type stars and found excess around 12% of them in the long-wavelength IRS band (30- 
34 /tm), but only 1% of the stars have detectable excess in the short wavelength band (8.5-12 /im). 

Exozodiacal clouds have been found around at least 20% of AFGKM nearby stars searched with 
IR interferometry (CHARA/FLUOR and Keck/KIN) (Absil et al., in prep.). There appears to be a weak 
positive correlation between the incidences of outer disks and exozodis (O. Absil, pers. comm). 



5.5 Cliemical Composition of Dust 

Mid- and far-IR spectra o f some debris disks reveal distinctive features dJura et all 120041: IChen et all 
I2OO6I iLawler et all l2009l) . which allows one to get insight into the mineralogy of the dust grains. 
For example, spectra of several disks were matched by a rnixture of amorphous and crystalline sili- 
cates, silica, and several other species dSchutz et all l2005l IBeichman et all l2005t iLisse et all l2007l 
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iLisse et al l l2008HLisse etalll2009l) . including possibly w ater ice (IChenetalll2008h . The most detailed 
study was perform ed for the systems HP 69830 ( Beichma n""etail 120051: iLisse et al.Ll2007l) . HD 11 3766 



jLisse et al. I l2008l) . and HD 172555 dLisse et all 120091) . However, interpretation of spectra is difficult 



and involves degeneracies, since the same spectra can be fitted with different mixtures of materials. 
Furthermore, laboratory spectra of various samples used as a reference have been obtained under con- 
ditions (grain sizes and temperature) that do not necessarily match those in real disks. Furthermore, the 
laboratory spectra ca n differ from one mea surement technique to another (e.g., KBr pellets vs. free- 
floating particles, see iTamanai et al.l l2006h . Despite these caveats, the results are useful. Beside the 
chemical composition, some constraints can be placed on dust grain sizes and morphology. For instance, 
the very presence of distinctive spectral features is indicative of small grains. 



6 DEBRIS DISKS: THINKING OF PLANETESIMALS 
6.1 Modeling "from the Sources" 

Results outlined in the previous section describe the "visible" dust component in debris disks, but they 
do not tell us directly about the dust parent bodies and dust production mechanisms. Therefore, many 
of the basic parameters of the debris disk remain obscure. For instance, while the cross-section area 
of the disk and thus the observed luminosity are dominated by small particles at dust sizes, the bulk 
of a debris disk's mass is hidden in invisible parent bodies and cannot be directly constrained from 
analysis of the dust emission. Equally, it remains unknown where exactly the planetesimals are located, 
although one expects that they orbit the star roughly where most of the dust is seen. Many properties 
of the planetesimals, such as their dynamical excitation, size distribution, mechanical strength, porosity 
etc. remain completely unclear. 

There is no direct way to infer the properties of invisible planetesimal populations from the observed 
dust emission. Dust and planetesimals can only be linked through models. This is done in two steps. 
First, collisional models can be used to predict, for a given planetesimal family (mass, location, age, 
etc.), the distribution of dust. Such models, both analytical and in the form of numerical codes, have 
become available in recent years (e.g. iThebault et al.L 120031: iKrivov et al. l l2006HThebauh & Augereaul 



120071: IWvatt eTan.l2007aLlLohne et al.ll2008l) . After that, thermal emission models, described in Sect.|5] 
have to be employed to compute the resulting dust emission. Comparison of that emission to the one 
actually observed would then reveal probable properties of the dust-producing planetesimal families. 

6.2 Constraints from Statistics of Debris Disks 

One way to constrain planetesimal properties with this method is to develop models for a long-term 
evolution of disks and to compare their predictions with the statistics of debris disks of different ages. 
As noted above, there is a general tendency of disk dustiness to decay with time, and this is attributed 
to the collisional depleti on of planetesimal famili es. In an attempt to gain theoretical understanding of 
the observed evolution, iDominik & DecinI (|2003|) assumed that dust is produced in collisions among 
equally-sized "comets". If this dust is removed by collisions, the steady-state amount of dust in such 
a system is proportional to the number of comets. This results in an M/Mq r/t dependence for 
the amount of dust and for the number of comets or the total mass of the disk. Under the assumption 
of a steady state, this result is valid even for more complex systems with continuous size distributions 
from planetesimals to dust. Tenuous disks, where the lifetime of dust grains is not limited by collisions 
but by transp ort processes like the Poynting-Robertson drag dArtvmowiczl Il997t iKrivov et all l2000l : 
IWvatll I2OO5I) . should follow M cx t^^ rather than M (x t~^. However, slo pes as steep as M oc 
are not consist ent with observat ions. Assuming a power- law dependence t~^. Greaves & Wvatd (l2003l) 



report ^ ^ O.S. lLiu et all (|2004 give 0.5 < C < 10, and lGreavesI (l2005b and lMoor et al.1 (|2006|) derive 



£ 1.0. Fits of the upper envelope of the distribution of luminosities over the age yield ^ « 1.0 as well 
jRiekeetaUl2005[) . This is an indication that the observed cold debris disks are collision-dominated. 
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I Wyatt et aP (l2007ah lifted the most severe simplifying assumption of the Dominik-Decin model, 
that of equal-sized parent bodies. A debris disk they consider is no longer a two-component system 
"comets + dust". Instead, it is a population of solids with a continuous size dis t ributio n, from planetes- 
imals down to dust. A key param e ter of t he description by iDominik & DecinI (l2003h is the collisional 
lifetime of comets, t. IWvatt et al] ( l2007al) replaced it with the lifetime of the largest planetesimals and 
worked out the dependencies on this parameter in great detail. Since the collisional timescale is inversely 
proportional to the amount of material, t cx 1/Mq, the asymptotic disk mass becomes independent of 
its initial mass (Fig. |5] top left). Only dynamical quantities, i.e. the disk's radial position and extent, 
the orbiting objects' eccentricities and inclinations, and material properties, i.e. the critical specific en- 
ergy and the disruption threshold, as well as the type of the central star determine the very-long-term 
evolution (Fig.|5] thin lines). 
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Fig. 5 Fractional lum inosity of dust arou nd a solar-like star as a function of age. Thick 
fines: analytic model of iLohne et"ai] (l2008b : thin lines: Eq. (20) of lWvatt etal] (l2007ah with 
= 3 X 10* erg (constant in their model), = 60 km,. Different panels demonstrate 
dependence on different parameters: ^/disk (top left), r (top right), dr/r (bottom left), and e 
(bottom right). A standard case with A/* = = 1, A/disk = lOMm, r = 30 AU , dr/r = 
1/2, and e = 0.10 is shown with solid lines in all panels. Adapted from lLohne et al\tOO^ . 



Most recently. iLohne etal] ( |2008[) lifted a major assumption in lWvatt et alj ( l2007al) that the critical 
specific energy needed for disruption is constant across the full range of sizes, from dust to the largest 
planetesimals. As a result, the size distribution of material is no longer a single power law. One con- 
sequence is that at actual ages of debris disks between ^10 Myr and ~'10 Gyr, the decay of the dust 
mass and the total disk mass follow different laws. The reason is that, in all conceivable debris disks, the 
largest planetesimals have longer collisional lifetimes than the system's age, and therefore did not have 
enough time to reach collisional equilibrium. The mass of visible dust at any instant of time is deter- 
mined by planetesimals of intermediate size, whose collisional lifetime is comparable with the current 
age of the system, with that "transitional" size gradually increasing with time. Under standard assump- 
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tions, the dust mass, fractional luminosity, and thermal fluxes all decrease as with ^ = 0.3 . . . 0.4 
(cf. Eq. (|9]l). However, specific decay laws of the total disk mass and the dust mass, including the value 
of ^, largely depend on planetesimal properties, such as the critical fragmentation energy as a function of 
size, the "primordial" size distribution of largest planetesimals, as well as the characteristic eccentricity 
and inclination of their orbits. 



6.3 Systems with Abnormally High Fractional Luminosity 

Bofli lWvatt" et al.l (12007 al) and lLohne et al.l (l2008l) models reproduce reasonably well the observed statis- 
tics of "cold" and many "warm" debris disks obtained with the Spitzer/MIPS instrument. However, there 
are individual systems that are not consistent with the models. One implication of the above models is 
that they predict a certain maximum fractional luminosity /max for a given age that could be produced 
by a debris disk collisionally evolving in a steady-state regime. Assuming physically reasonable initial 
disk masses of ^ SOA/^ (corresponding to the mass of solids in the minimum-mass solar nebula), /max 
of debris disks around solar-type stars older than ^ 1 Gyr would be as low as ^ 10~^ (Fig. |5] top 
left). But some known debris disks have / ^ /max- These are "warm " disks within ~ 10 AU, exem- 
plified by the de bris disk of an A star L eo (iMoerchen et"an. l2007[) . as well as some disks o f FGK 
stars: HP 2 3514 (Rhee et al.', '20081), rj Crv dWvatt et al.ll2005l) . and HD 69830 (iBeichman et al.ll2005l 



iLisse et al.l 12007; Pa yne et a l., 2001). The same problem arises for exozodiacal clouds, where /max is 
quite low because the collisional evolution as close as at ^ 1 AU from the sta r must be very rapid. Thus 
it is not a surprise that for instance, the exozodi of Vega (lAbsil et al.l l2006l) appears "too dusty" to be 
explained with a steady-state collisional cascade in an a steroid belt analog. The origin of the excessive 



warm dust and exozodis is currently a matter of debate dWvatt et"ai] , l2007al IWvattl 20081 iPavne et al 



120091 IWvattetaD. 120101) . The obse rved dust may be debris from recent collisions between two large 
planetesimals or planetary embryos (iRhee et al.ll2008l) . Al ternatively, it can be transported inward from 
the outer Kuiper belt through one or another mechanism dWvatt et al.l l2007al) . It is possible that some 
of the systems with a high amount of warm d ust currently undergo short-lasting, LHB-like instability 
phase jWvatt et al.L l2007al iBooth et all |2009|) . Finally, systems like rj Crv may formally be explained 
with a collisional cascade in a planetesimal swarm with extremely ecc entric orbits (e ~ 0.99), although 
the origin of such population remains challenging dWvatt et al.ll2010l) . 



6.4 Probing Planetesimal Formation Mechanisms 



Since the formation of debris disks is a natural by-product of planetesimal and planet accumulation 
processes, debris disks can help distinguishing between possible planetesimal formation mechanisms. 
The classical scenario is based on a more or less continuous growth of objects from micron-sized grains 
through bigger aggregates that decouple from the gas to, potentially, pla netary embryos. Especially 
the last stage has been studied to predict planetesimal distributions (e.g. iKenvon & Bromlevl l2008l 
and previous work). However, theoretical considerations have revealed obstacles such as the so-called 
meter barrier. The latter is a "double bottle- neck" of the grow'th: first, meter-sized obje cts should be 
lost to the central star as a result of gas drag dWeidenschillind, [l977t iBrau er et al.l ^2007*). and second, 
further agglomeration of meter-sized objects upon collision is problematic CBlum & Wurm. 200^. So, 
in the last years, a competing scenario was invoked that circumvents these barriers and forms large 
planetesimals directly through local gravitational instability of sol ids in turbulent gas disks, which can 
be triggered, for instance, by trans ient zones of high pressure (Jo hansen et al.l l2006b or by streaming 
instabilities (jjohansen et al.l l2007b . Both scenarios result in distinctly different initial conditions for 
debris disk evolution: the former predicts a substantial population of smaller km-sized planetesimals, 
while the latter leads to a top-heavy size distribution dominated by objects of roughly 100 km in size. 
As both scenarios would result in different dustiness of a debris disk at a given age, there appears to be 
a principal possibility to distinguish between them by confronting the model predictions to debris disk 
observations (Lohne et al., in prep.). 
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Constraining the top end of the planetesimal size distrib ution is also importan t to understand the 
mecha nism that stirs debris disks. Extending the formalism of lWvatt et all (l2007albl) . [Kennedy & Wv"at3 
(l2009l) devised an analytic model for stirring of debris disks by the formation of Pluto-sized objects. 
They showed that the 24 and 70 i-im statistics of disks around A-type stars is consistent with this self- 
stirring scenario. That also demonstrated that, to satisfy the constraints from observations such as the 
early "rise and fall" of peak 24 /im excesses around ^ 10 Myr (Sect.|5]), the planetesimal belts must be 
reasonably narrow (dr/r ^ 0.5) and have a minimum radius of ~ 15 AU. 

6.5 Constraints from Debris Disks with Known SEDs 

As we could see, useful constraints on the planetesimal properties can be placed for the statistics of 
many debris disks, even though observational data on each of them can be limited to a few photometric 
points, such as 24 and 70 yum Spitzer/MIPS measurements. For a limited number of disks, for which 
photometry points sufficiently pro be the entire SED f rom mid-IR to sub-mm or mm wavelengths, ad- 
ditional constraints can be placed. [Rrivov et all ( |2008|) analyzed five G2V stars with good data (IRAS, 
ISO/ISOPHOT, Spitzer/IRAC, /IRS, /MIPS, Keck II/LWS, and JCMT/SCUBA). For all five systems, 
the data points could be reproduced within the error bars with a linear combination of two planetesimal 
belts (an "asteroid belt" at several AU and an outer "Kuiperbelt"). This automatically gave the desired 
estimates of planetesimals (location, total mass etc.). In particular, the cold emission (with a maximum 
at the far-IR) is compatible with "large Kuiper belts", with masses in the range 3-50 Earth masses 
and radii of 100-200 AU. These large sizes trace back to the facts that the collisional model predicts 
the observed emission to stem from micron-sized dust grains, whose temperatures are well in excess 
of a b lackbody temperature at a given distance from the star (as discussed, e.g., in [Hillenbrand et all 
[20081 see also Fig. I?]). This conclusion is rather robust against variation in parameters of collisional and 
thermal emission models, and is rou ghly consistent with disk radii revealed in scattered light images 
(e.g. HD 107146. [Ardila et all [20041) . Still, quantitative conclusions about the mass and location of the 
planetesimal belts would significantly depend on (i) the adopted model of collision outcomes (which, 
in turn, depend on the dynamical excitation of the belts, i.e. on orbital eccentricities and inclinations 
of planetesimals) and (ii) the assumed grains' absorption and emission efficiencies. For example, a less 
efficient cratering (retaining more grains with radii ~ 10 /im in the disk, see Fig. [2] left) and/or more 
"transparent" materials (making dust grains of the same sizes at the same locations colder) would result 
in "shifting" the parent belts closer to the star. 

6.6 Constraints from Resolved Debris Disks 

The best constraints of planetesimals and their collisio nal evolution can b e found from a combined 
analysis of the SED and resolved images. Most recentlv. [Miiller et al.l ([2010l) analyzed the archetypical 
debris disk of Vega. It has been argued before that the resulting photometric data an d images may b e 
in contra diction with a standard, steady-state collisional scenario of the disk evolution iSu et al.l [2005f). 
However, [Miiller et al.[ ( [2010l) were able to reproduce the spectral energy distribution in the entire wave- 
length range from the near-IR to millimeter, as well as the mid-IR and sub-mm radial brightness profiles 
of the Vega disk with this scenario (Fig.|6]l. Thus their results suggest that the Vega disk observations are 
not in contradiction with a steady-state collisional dust production, and they put important constraints 
on the disk parameters and physical processes that sustain it. The total disk mass in ^ 100 km-sized 
bodies is estimated to be ~ 10 Earth masses. Provided that collisional cascade has been operating over 
much of the Vega age of ^ 350 Myr, the disk must have lost a few Earth masses of solids during that 
time. It has also been demonstrated that using an intermediate luminosity of the star between the pole 
and the equator, as derived from its fast rotation, is required to reproduce the debris disk observations. 
Finally, it has been shown that including cratering collisions into the model is mandatory. 

The appearance of the Vega disk in resolved images, with a ring of large particles seen in sub-mm 
and extended sheet of small grains observed in the far- and mid-IR, conforms to the standard debri s 
disk model described in Sect. [D Some other resolved disks, such as that of HR8799 dSu etali[2009l) . 
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Fig. 6 Two models of the Vega disk jMiiller et al.l I2OIOI) and their comparison with observa- 
tions. The "reference model" assumes = 28L0 and a ring of planetesimal with semimajor 
axes from 80 to 120 AU and maximum eccentricities of 0.1. The "best-fit model" adopts 
Li, = 45^0, a planetesimal ring from 62 to 120 AU with maximum eccentricities of 0.05, 
and a slightly steeper distribution of fragme nts in binary collisions. The reason to vary th e 
stellar luminosity is that Vega is a fast rotator dPeterson et al.l I2OO6I : lAufdenberg et"alll2006l) . 
which makes stellar parameters functions of the stellar latitude. Top left: grain size dust distri- 
bution (cross section area density per size decade) in the "reference modefXthick solid line) 
and the "best-fit model (dashed line), both in the center of the initial planetesimal ring. Note 
our using 2tis instead of s: since particles with the size parameter 2tis/\ ^ 1 emit most effi- 
ciently, 27rs roughly gives the typical wavelength of the emission. This alleviates comparison 
between the size distribution and the SED (bottom left). Top right: radial profile of the normal 
optical depth for the same disk models. Bottom left: corresponding SED in the form of the 
excess ratio (ratio of dust flux to photospheric flux). Symbols with error bars are data points 
(large square, large circle, and large triangle mark 24 /ini, 70 /im, and 160 /im excess ratios 
deduced from Spitzer/MIPS images). Bottom right: modeled (lines) and observed (symbols, 
same as in the bottom left panel) surface brightness profiles at 24, 70, and 160 /^m. The shaded 
areas around the data points correspond to their error bars. Adapted from \Muller et al.\ iZOld) . 



look similar However, the disks of two other A-type stars, Fomalhaut dStapelfeldt et al.l l2004 and 
/3 Leo (Stock et al., in prep.), appear narrow at Spitzer/MIPS wavelengths. The existence of two types 
of disks - wide ones with fuzzy sensiti vity-limited outer edges and narrow ones with w ell-defined outer 
boundaries — was also concluded by iKalas et all (|2006|) from scattered-light images. iThebault & Wul 
d^OS) showed that disks with sharp outer edges may arise from a parent belt of planetesimals with very 
low eccentricities (e ~ 0.01 or less). However, this would imply the absence of planetesimals larger than 
several tens of meters, as these would stir the parent belt to higher eccentricity values. It is challenging 
to explain why, for instance, the protoplanetary disk of Fomalhaut that has built at least one massive 
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planet at - 100 AU dKalas et all I2008I) . but failed to produce at least kilometer-sized planetesimals 
only slightly farther out from the star. 



7 DEBRIS DISKS: THINKING OF PLANETS 

7.1 Structures of Debris Disks as Indicators of Planets 

The usefulness of debris disks goes beyond probing the history of planetesimal formation and planetes- 
imal properties: they can also point to planets. Indeed, debris disks gravitationally interact with - often 
suspected, but as yet undetected - planets at every moment in their life, and various footprints of those 
interactions can be evident in observations of debris dust. As mentioned above, planets may stir plan- 
etesimal disks, launching the collisional cascade. Stirring by planets may be comparable in efficiency 
with stirring by largest embedded planetesimals of roughly the Pluto size. Even close-in plan ets can 
suffic i ently stir the outer disk , if their orbits are eccentric and the systems are sufficiently old dWvatli 
I2OO5I : iMustill & Wvatli I2OO9I). For instance, the kno wn radial velocity planet of e Eridani could have 
ignited its "cold" debris disk (IMustill &Wvattl 120091) . 

Besides the stirring, a variety of more direct "fingerprints" of planets can exis t in the disks. Inner 
gaps tens of AU in radius commonly observed in resolved debris disk im ages (e.g. Moro-Martm et all 



I2OO7I and references therein) and inferred from the debris disk statistics dKennedv & Wvatf 20091) ^ 



likely the result of clearance by planets dPaber & Ouillenl |2007|) (although they ma y altern atively be 



ascribed to preferential planetesimal formation at larger distances in some scenarios. iRice et a l. 2006) 



Accordingly, dedicated quests of planets in systems with known "punched" debris di sks by using direc t 
imaging from space or with adaptive optics have been initiated in the recent years dApai et all l2008l) . 
The direct imaging method is technically challenging, and only has chances for success in relatively 
young systems with ages younger than a few hundred of Myr, where the planets have not yet cooled 
down, but even in that case current searches are not yet sensitive to planets below about one J upiter mass. 
Desp ite that, first pl anets in the inner vo ids have recently been found around FomaUiaut dKalas et all 
I2OO8I) and HR 8799 dMarois et alll2008l) . 

Another indication for planets could be a clumpy disk structure observed, for instance, around Vega 
dWvattLl2003l ) and e Eri dOuillen & Thorndikell2002l). Clumps may arise from plane tesimal capture into 



mean motion resonances with unseen planets (e.g. Wvatil2006llKrivov et alll2007l) or, in systems with 



a str ong P-R or stellar wind drag, from capture of dust grains drifting inward into such resonances 
(e.g. lLiou&ZookLll999t lM oro-Martm & Malhotrail2002l) . A cave at is that individual brea kups of large 
planetesimals could sometimes provide an alternative explanation dGrigorieva et all [20071) . Attempts to 
find planets in the systems of Vega and e Eri by direct imaging have not yet been successful ditoh et all 
I2OO6: .Mare ngo et al, 2006; Janson et al, 2007, 2008). 

Warps, such as those observed in the /3 Pic disk, could be due to secular perturbations by a planet in 
an inclined orbit dMouillet et allll997l:lAugereau et all 1200 lb . In fact, there are many more phenomena 
in (3 Pic disk that all point out to at least one Jupiter-mass planet at ^ 10 AU (see iFreistetter et all 
|2007L for a summary). T he predicted /? Pic planet may have already been detected with VUT/NACO 



by Lagrange et alj d2009l). alth ough this result still remains to be confirmed with a second-epoch de- 
tection dFitzgerald et al 1 120091) . However, large-scale asymmetries o bserved in resolved debris disks. 



especially the wing asy mmetries in the disks of /? Pic ( iMou illet et an,ll997l). HD 32297 dKalasl I2OO5I: 



iFitzgerald etall l2007l) . and HD 15115 ("the blue needle". iKalas et al.'. 2007^ can be attr ibuted to the 
disk i nteraction with the ambient interstellar dust dArtvmowicz & Clampini, 1 19971) or gas dOebes et all 
l2009h . 
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7.2 Incidences of Debris Disks and Planets around the Same Stars 



Current observational statistics shows that at least 7% of stars searched with radial velocity (RV) method 
have planets ranging from super- Jupiters down to super-Earths (Udry & Santos, 2007J3- On the other 
hand, IR excesses over the photospheric flux indicative of (cold) debris di sks have been found for about 
15% of main-sequence (MS) stars with spectral c l asses from A to K (|Su et"ail l2006t ISiegler et all 
120071; iTrilUng et all l2008t [Hillenbrand eTall l2008h . lOreaves et all (I2OO4I) were apparentlv the first to 
raise natural questions: "do stars with planets also have disks?" and conversely, "do stars with disks 
also have planets?". Answering these questions is hampered by the fact that the incidences of both 
planets and disks, retrieved from observational surveys, are biased in various ways, depend on the cur- 
rent precision or sensitivity of instruments used, and thus may not reflect the reality. It is not a sur- 
prise, therefor e, that the stati s tical interrelation between the presence of planets and cold disks remains 



controversial ( Greaves et all 20041; Beichman et all 

■ ml 



;n the presence or planets and cold disKs remams 
20051; iMoro-Martin et all 12007 al: lOreaves et all 



l2006l:lBeichman et all l2006t iKosp^l et alll2009HBrvden et alll2009al) . Neverthel ess, it has been found 



that the incidences of debris disks do not correlate with the stellar metallicity dGreaves et all 120061 ; 
iBeichman et"!!! l2005t ISrvden et alll2009ah . wh ereas the typical me tallicity of stars with known plan- 
ets is well known to be greater than solar (e.g. I Wright et all |2004|) . This might suggest the absence 
of correlation between the presence of planets and presence of debris disks, but this conclusion may be 
premature. For instance, the incide nces of low-mass pl anets ("Neptunes"), similar to debris disks, do not 
seem to correlate with metallicity dSousa et all l2008h . And it is not known if such a correlation exists 
for more distant planets that cannot be found with the RV techniques. Thus a correlation between the 
presence of planets and presence of debris disks cannot be ruled out. In fact, there are hints for the de- 
bris disk harboring systems with a planet detection to poss ess brighter, or dustier, disks thus higher disk 
detection rates) than systems without a planet detection dBeichman et al.l 120051: iTrilling et all l2008l 
iBrvden et alll2009ah . 



8 DEBRIS DISKS: THINKING OF PLANETARY SYSTEMS 

We now focus on those planetary systems that, like our solar system, are known to contain, apart the 

central star and one or more planets, also an outer and/or inner planetesimal belt, or a combination of 

them. A handful of "full" systems in this sense have been foun d so far, but only fir st steps have been taken 

to elu cidate th eir detailed s t ructur e. Exampl es are Fomalhaut (Kalas et alll2008l) . e Eri dBackman et al 
I 1 ^ . I, . . ( I I - ^- rr^T"!^ iTT^^^^^^^^^^^^^^^ 



20091) 



Eri dLiseau et al .','2008'). HR8799^^Reidemeister et al.'.'2009'). HP 38529 ([Moro-Martm et al 



2007bl) . and HP 69830 dBeichman et all 12005; Lisse et al, 2007; Payne et al... ,2009 

As an example, we take the system HR 8799, a nearby A-type star with three directly imaged 
planetary candidates (Marois et al.', '2008*), of which the outermost one was recently confirmed by 
[Lafreniere et al. (2009) and by Fukagawa et al. (2009) by analyses of the archival data taken earlier 
Apart from the planets, HR 8799 has long been kno wn to harbor cold circumstellar dust responsible 
for excess emission in t he far-IR discovered by IRAS dSadakane & Nishidaill986tlZuckerman & Songl 
120041; iRhee et alll2007|) The r ather strong IR excess has also been confirmed with ISO/ISOPHOT mea- 
surements ( Moor et all 20061) , and recently the dust disk was resolved with Spitzer/MIPS at 24 and 



70 /im by ISu et al. Additionally, Spitzer/IRS measurements provided evidence of warm dust 



emission in the mid-IR dluraetaU 120041 IChen etalll2006h. Both cold and warm dust emission is in- 
dicative of two dust-producing planetesimal belts. Reidemeister et al] d2009b undertook a coherent anal- 
ysis of various observational data for all known components of the system, including the central star, 
imaged companions, and dust. This has led to a view of a complex planetary system with an asteroid- 
like dust producing inner belt at ~ 10 AU, three planets of few Jupiter masses likely locked in a triple 
Laplace resonance, an outer Kuiper-like planetesimal belt at ~ 100 AU, associated with a tenuous dust 
disk that extends to several hundreds of AU from the star. Although many parameters of this system are 



' True fractions ai'e likely mucli higher. Most recent analysis of the HARPS GTO survey data suggests that 30 ± 10% of solar- 
type stars may host close-in (periods < 100 days), low-mass (< SOAfg) planets (Udry et al., paper presented at the ESO/CAUP 
conference "Towards Other Earths: perspectives and limitations in the ELT era", Poito, Portugal, 19-23 October 2009). 
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very different from those of our solar system (young age, luminous central star, very massive planets), 
the arrangement of these components appears similar (Fig.|7h)- 

a) HR 8799 b) HD 69830 




lOOAU 50AU 



Fig. 7 Schematics of four planetary systems with at least one planet and one planetesimal 
belt: (a) HR 8799, (b) HD 69830, (c) e Eri, and (d) HD 38529. Only known — not just pre- 
sumed — components are shown. Sketches are not to scale. Distance rulers roughly reflect the 
size of the entire systems, the inner parts of the e Eri and HD 38529 are artificially enlarged. 
The planets are marked with filled circles whose size is representative of their masses. 



Other "full" planetary systems show less similarity to the solar system. For instance, the system of 
HD 69830, a nearby 4-1 Gyr old KOV star , has a very compact configuration with three Neptune-mass 
planets within ~ 1 AU (iLovis et al.l 2006|) an d a "warm" (excess at 8-35 /^im) dust belt just exterior 
of the planetary region jBeichman et al. , 20051) (Fig.|7]5). A natural guess that dust can be collisionally 
produced in an asteroid belt analog fails, as the dust fractional luminosity is too high (see Sect. |6]l. 
Several scenarios have been pr oposed to explain th e origin of the warm emission, but all of them have 
so far unresolved problems (see lPavne et al 1 120091 for a detailed discussion). 

The nearby (3.2 pc ), ^ 1 Gyr old K2 V star e E ridani has a ring of cold dust at ~ 60 AU seen in re- 
solved sub-mm images dGreaves et all 19981 2005b. w hich is encompassed by an extended disk of warm 
dust resolved by Spitzer/MIPS ( Backman et all 20091) . The star is orbited by an RV planet dHatzes et all 
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\200Ci iBenedict et all 120061: iButler et all l2006h with a semimajor axis of 3.4 AU, and another outer 
plan et is expected t o orbit at ^ 40 AU, producing the inner cavity and a clumpy structure in the outer 
ring dLiou & Zooklll999l:IOzernoy et al.ll2000HOuillen & Thorndike ', 2002; Deller & Maddisonll2005b 
(Fig Et)- The excess emission at A ^ 15 /xm in a Spitzer/IRS spectrum (Backman et al.„i2009l) indicates 
that there is warm dust at just a few AU from the star Its origin is obscure , as an inner "aste roid belt" 
that could produce this dust would be destroyed by the known inner planet ( Brogi et"aill2009l ). 

One m ore example of a " f ull" pla netary system is the one around a 3.5 Gyr old G8III/IV star 
HD 38529 jMoro-Martm et all l2007bt) . It hosts two known RV companions with M sini of 0.8 and 
12.2 Jupiter masses, semimajor axes of 0.13 AU and 3.74 AU, and eccentricities of 0.25 and 0.35, 
respectively. Spitzer observations show that HD 38529 has an excess at 70 /im, which is compatible 
with dust coming from a "Kuiper belt" between 20-50 AU (Fig. |2}1). 

We do not know yet which combination of planets and dusty planetesimal belts is common, how 
these components are typically arranged, and which circumstances are decisive to set one or another 
type of architecture. However, it is particularly interesting to see how the dust helps to better understand 
the systems. We take an example of HR 8799 again. Fi rst, dust helped to narrow the previou sly reported 
range of age estimates of the system, from 30Myr ( iRhee et al.l l2007h to w 1100 Myr (ISong et al.l 
l200lb . Namely, the measure d IR excess rati o of ^ 100 at 60 -90 fim would be typical of a debris disk 
star of age < 50 Myr (see ISu et al. 1 120061 their Fig. 5). Good age estimates are crucial from every 
point of view, but they also improved the mass estimates of three planets. Indeed, the masses of three 
planets are estimated from their observed luminosities with the aid of the evolutionary models, and the 
results strongly depend on the age of the system. Old ages would inevitably imply companion masses 
in the brown dwarf range, and three co mpanions would then be dynamically unstable, but for the ages 
of ^ 50 Myr there are stable sol utions (iGozdziewski & Migaszewskl l2009t iReidemeister et all l2009l 
iFabrvckv & Murrav-Clavl l2010l) . Second, the fact that dust in the outer ring was detected as close as 
at 120 AU from the star and that dust in the inner ring extends to at least 10 AU away from it, sets 
independent upper limits of planetary masses that are consistent with masses from evolutionary models. 
Third, dust helped to constrain the orientation of the entire system, whose knowledge is needed to 
convert the sky-projected positions into their true position is space, setting correct initial conditions 
for the dynamical stability analyses. IReidemeister et all (l2009b showed that stellar rotational velocity is 
consistent with an inclination of 13 - 30°, whereas i ^ 20° is needed for the system to be dynamically 
stable, thus arguing for a probable inclination of 20 - 30°. T hat nearly face-on orientation was confirmed 
by resolved images of the outer dust disk obtained bv lSu et al.. (.2009.) . 



9 OPEN QUESTIONS 

After twenty-five year s of research on circumstellar deb r is disk s since the first discoveries around Vega 
(lAumann et al.l 1 19841) and /3 Pictoris dSmith & Terrilel 1 19841) astronomers can unambiguously place 
them into the general context of formation and evolution of planetary systems. We know that stars 
form in collapsing interstellar clouds, that young stars are surrounded by dense, gas- and dust-rich 
protoplanetary disks in which planets can form, and that debris disks is what remains around the star 
after the dispersal of that primordial gas at ages of ~ 10 Myr or slightly below. This view is directly 
supported by the observations that show a clear difference between the disks around younger and older 
stars. Below those ages, the mass of dust in circumstellar disks retrieved from the IR excesses is typically 
~ lO^Afg (Earth mass), and above them it drops to below one ~ IM^. Applying a standard gas-to-dust 
ratio of ~ 100:1 to protoplanetary disks, their mass is estimated to be ~ 10^ M^, which is consistent 
with direct gas detections. Direct detection of gas in debris disks is very difficult, and determination of 
its mass rather uncertain even where it was detected, but is it likely that the gas mass does not exceed a 
fraction of the Earth mass even in the gas-richest debris disks such as /3 Pic and HD 32297. 

One concludes that a certain star possesses a debris disk if thermal emission of circumstellar dust 
and/or stellar light scattered by it has been detected, and if the am ount of dust is below a certain limit, 
which is taken to be 10^^ in terms of the fractional luminosity dUagrange et al.l l2000l) . The latter is 
required to distinguish between debris disks and protoplanetary disks. An additional criterion is that the 
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central star should be older than several Myr. From observations, we can only say that dust is present, 
and from resolved images we know that it is present in the form of a disk rather than an envelope, 
but a question is where that dust comes from. It is easy to estimate the dust lifetimes around a star 
against collisions, P-R or stellar wind drag, or destructive processes such as photosputtering. Which 
removal mechanism is dominant depends on the properties of the central star, dustiness of the disk, and 
dust composition — but the resulting dust lifetime is usually well below 1 Myr. Since debris disks are 
commonly observed around stars with ages up to several Gyr, the probability that we are "catching" 
freshly-born dust is vanishingly low, and so, one or another mechanism to continually replenish the 
disks is required. There is currently little doubt that the main mechanism that supplies dust is collisions 
between planetesimals that orbit the debris disk host star Thus the very existence of debris disks is 
a strong evidence for planetesimal populations orbiting a significant fraction of main sequence stars. 
This, in turn, tells us that planetesimal formation in protoplanetary disks is efficient. However, we do 
not know yet how these planetesimals typically form. Future debris disk studies may help discerning 
between alternative mechanisms of planetesimal formation, since they predict different size distributions 
of planetesimals and/or preferential formation of planetesimals in different regions of protoplanetary 
disks, which may be reflected in dust emission at the later, debris disk evolutionary phase. 

What is the true incidence of debris disks around main-sequence stars? Frequency of dete ction of 
protoplanetary disks in clusters younger than ^ 1 Myr is almost 100% dHernandez et aklboOTl) . but the 
detection fraction of debris disks is w 15% across the spectral stars and ages (but « 60% for A-stars 
at ^ 10 Myr). However, this may only reflect the sensitivity of current observations. It is possible that 
many other stars, perhaps even 100%, harbor fainter debris disks that as yet remain undetectable. The 
same applies to planets. Current RV detections lead to frequencies of Jupiter-mass planets of ^ 7%, 
and of lower-mass planets of ~ 30%, but many more — perhaps nearly all — stars may harbor planets 
that simply escape detection with present-day techniques and instruments. Statistical relations between 
systems with planets and disks remain unclear We only know that there are stars that have both. Thus, 
the typical questions are: Do nearly all mature stars harbor planets, planetesimals, and dust? Are there 
systems that have planets, but do not have planetesimal belts? Are there systems with planetesimal belts 
but without planets? And are there systems, in which protoplanetary disks failed to produce both? 

Accordingly, it is not yet possible to judge about the place of our own solar system's planetesimal 
famihes, including transneptunian objects, asteroids, and comets, in the menagerie of debris disks in 
other planetary systems. It is clear that the solar system's dust disk, dominated by dust in the Kuiper belt 
region, is by orders of magnitude fainter than known debris disks of other stars. But it is not yet clear 
why, and, due to current observational limitations, it is not yet known if any other stars possess disks as 
tenuous as ours. And, if they do, what is more common — brighter or fainter disks? 

We saw that statistics of IR excesses of hundreds of debris disk stars and their evolution with sys- 
tem's age, as well as observed SEDs and resolved images of some brightest debris disks appear con- 
sistent with a standard collisional evolution scenario. It implies one (or more) planetesimal belt with 
a steady-state cascade of successive collisions that grinds the planetesimals all the way down to dust, 
whose distribution is then shaped by gravity and radiation pressure forces. In such systems, steady- 
state collisional models place useful, albeit often weak, constraints on size and spatial distributions of 
planetesimals, their dynamical excitation, and mechanical strength. However, it remains unclear what 
mechanism stirs planetesimals to the extent that allows the collisional cascade to operate. Are these 
largest, Pluto-sized planetesimals embedded in the disks, neighboring planets, or both? On the other 
hand, there are systems — or their parts — that are certainly not consistent with a steady-state colli- 
sional evolution. These are exemplified by systems with substantial "hot" dust emission in the mid-IR, 
such as HD 69830 and rj Crv, and by systems with "exozodiacal" dust within ~ 1 AU from the star. 
Here, there is too much dust to be explained with dust produced locally in a steady-state collisional 
cascade. It remains unclear whether non-collisional mechanisms such as sublimation, one-time major 
collisions, or non-local production in outer belts with subsequent inward transport are at work in such 
systems. 

Further open questions are related to the disk structure. Most, if not all, known debris disks have 
inner cavities suggesting that planetesimals are not present there. Nonetheless, it is not known whether 
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the lack of planetesimals is because they originally fail to form there or were removed by planets later. 
Interpretation of azimuthal and vertical structures such as clumps, spirals, warps, and wing asymme- 
tries remains ambiguous, too. They can be naturally explained by gravitational interactions of solids 
with planets, but many can equally well be attributed to "non-planetary" mechanisms such as major 
coUisional break-ups, stellar flybys, or interaction with the surrounding interstellar gas and dust. 

We finally note that detailed studies of thermal emission of debris dust, notably the analyses of 
high-resolution mid- and far-IR spectra, offer a unique way to probe chemical composition of dust and 
its parent planetesimals. Only the first steps have been taken in this direction. For instance it is not yet 
clear, why some of the disks reveal spectral features and some others do not, to which extent the inferred 
composition of dust reflect that of its parent planetesimals, and whether the dissimilarity of the observed 
spectra traces back to different chemical evolution and mixing scenarios at the preceding, protoplanetary 
disk stage. Additional information can also be gained by analysis of gas in young debris disks, which is 
thought to be a product of dust grain disintegration by various mechanisms. 

Future observations and modeling work should bring answers to these and other questions. On any 
account, debris disks must be treated as an important constituent part of planetary systems, and there is 
little doubt that they can serve as tracers of planetesimals and planets and shed light on the planetesimal 
and planet formation processes. Thus the lesson for planetary system research is: "follow the dust". 
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